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Abstract 

The energy of a simple graph G, denoted by E{G), is defined as the sum 
of the absolute values of all eigenvalues of its adjacency matrix. Denote 
by Cn the cycle, and the unicyclic graph obtained by connecting 
a vertex of Cq with a leaf of Pn-e ■ Caporossi et al. conjecture that 
the unicyclic graph with maximal energy is P^ for n = 8, 12, 14 and 
n > 16. In"Y. Hon, I. Gutman and C. Woo, Unicyclic graphs with 
maximal energy. Linear Algebra Appl. 356(2002), 27-36", the authors 
proved that E{P^) is maximal within the class of the unicyclic bipartite 
n- vertex graphs differing from Cn ■ And they also claimed that the 
energy of Cn and is quasi-order incomparable and left this as an 
open problem. In this paper, by utilizing the Coulson integral formula 
and some knowledge of real analysis, especially by employing certain 
combinatorial techniques, we show that the energy of Pn is greater than 
that of Cn for n = 8, 12, 14 and n > 16, which completely solves this 
open problem and partially solves the above conjecture. 
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1 Introduction 



Let G be a simple graph of order n, A{G) the adjacency matrix of G. The characteristic 
polynomial of A{G) is usually called the characteristic polynomial of G, denoted by 

x) = det(x/ - A{G)) = + aix""^ H h a„, 

It is well-known [2] that the characteristic polynomial of a bipartite graph G takes the 
form 

Ln/2J ln/2\ 
k=0 k=0 
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where = (— l)'^a2fc and 62^ > for all fc = 1,..., Yn/2\, especially = ao = 1. 
Moreover, the characteristic polynomial of a tree T can be expressed as 

\n/2\ 
fc=0 

where m(T, /c) is the number of /c-matchings of T. 

For a graph G, Let Ai, A2, . . . , A„ denote the eigenvalues of its characteristic polyno- 
mial. The energy of a graph G is defined as 



i=l 



This definition was proposed by Gutman [4]. The following formula is also well-known 

E{G) = -r ^log|xXG,Va;)|dx, 

where = —1. Furthermore, in the book of Gutman and Polansky [7], the above equality 
was converted into an explicit formula as follows: 



1 /'+°° 1 



' Ln/2J \ " / Ln/2J ^ 

fc=0 / \ fc=0 



dx. 



For more results about graph energy, we refer the reader to the recent survey of Gutman, 
Li and Zhang [6]. 

For two trees Ti and T2 of the same order, one can introduce a quasi order -< in the 
set of trees, namely, if m{Ti,k) < m(T2,k) holds for all k > 0, then define Ti ^ T2, 
and so Ti ^ T2 implies E{Ti) < E{T2) (e.g. [3]). Similarly, one can generalize the quasi 
order to the cases of bipartite graphs (e.g. [14]) and unicyclic graphs (e.g. [8]). The 
quasi order method is commonly used to compare the energies of two trees, bipartite 
graphs or unicyclic graphs. However, for general graphs, it is difficult to define such a 
quasi order. If, for two trees or bipartite graphs, the above quantities m(T, k) or |afc(G)| 
can not be compared uniformly, then the common comparing method is invalid, and this 
happened very occasionally. Recently, for these quasi-order incomparable problems, we 
find an efficient way to determine which one attains the extremal value of the energy, see 
[10-13]. 

Let Gn be the cycle, and be the unicyclic graph obtained by connecting a vertex 
of Gq with a leaf of P^-e • In [1], Caporossi et al. proposed a conjecture on the unicyclic 
graph with the maximum energy. 

Conjecture 1.1. Among all unicyclic graphs on n vertices, the cycle Gn has maximal 
energy if n < 7 and n = 9,10,11,13 and 15. For all other values of n , the unicyclic 
graph with maximal energy is ■ 

Theorem 1.2. Let G he any connected, unicyclic and bipartite graph on n vertices and 
G^Gn. Then E{G) < E{P^) . 



In [9], the authors proved Theorem 1.2 that is weaker than the above conjecture, 
namely that E{P^) is maximal within the class of the unicyclic bipartite n-vertex graphs 
differing from C„ . And they also claimed that the energy of Cn and is quasi-order 
incomparable. In this paper, we will employ the Coulson integral formula and some 
knowledge of analysis, especially by using certain combinatorial techniques, to show that 
E{Cn) < E{P^), and then completely determine that P^ is the only graph which attains 
the maximum value of the energy among all the unicyclic bipartite graphs, which partially 
solves the above conjecture. 

Theorem 1.3. For n = 8, 12, 14 and n > 16, E{P^) > E{Cn). 



2 Main results 

We recall some knowledge on real analysis, for which we refer to [15]. 
Lemma 2.1. For any real number X > —1, we have 

X 



l + X 



< log(l + X)<X. 



The following lemma is a well-known result due to Gutman [5], which will be used in 
the sequel. 

Lemma 2.2. If Gi and G2 are two graphs with the same number of vertices, then 



1 



"+00 



E{G^) - E{G2) = - log 



TT 



ix) 



(f){G2-,ix) 



dx. 



In the following, we list some basic properties of the characteristic polynomial (j){G, x) , 
which can be found in [2]. 

Lemma 2.3. Let uv be an edge of G . Then 

(p{G, x) = (j){G — uv, x) — (j){G — u — v,x) — 2 0(G — C, x) 

cec{uv) 

where C{uv) is the set of cycles containing uv . In particular, if uv is a pendent edge with 
pendent vertex v , then (l>{G, x) = x 4>{G — v,x) — (piG — u — v,x) . 

Now we can easily obtain the following lemma from Lemma 2.3. 
Lemma 2.4. 4>{Plx) = x(P{Pti,x)-4>{Pt2,x) and a;) = 0(P„, x)-</)(P„_2, a;)-2. 

By some easy calculations, we get 0(P|,x) = x* — 8x^ + 19x^ — 16x^+4 and (j){P^,x) = 
x^ — 7x^ + 13x'^ — 7x. Now for convenience, we introduce some notions as follows 



^ yi(x)(x^-l)-x ^ Y2ix)ix' - 1) - X 

^'^""^ {Y,{x)Y - Y,ix) ' "^'^"^^ iY2{x)y - Y2{x) ' 

yi(x)0(p|,x)-0(p6,x) y2(x)0(p|,x) -0(P7^x) 



iY,{x)f - (Fi(x))7 ' {Y2{x)f - {Y2{x))'' 



It is easy to verify that Yi{x) + Y2{x) = x, Yi{x)Y2{x) = 1, Yi{ix) = s+^^^IMi and 
Yo(ix) = We define 



f8 = x^ + 8x^ + 19a;^ + 16x^ + 4, /t = + 7x^ + 13x^ + 7x 



and 



Zi[x) = -lYiiix) = , Z2{x) = -iY2(ix) = — 



Lemma 2.5. For n > 10 anc? x ^ ±2, the characteristic polynomials of and Cn have 
the following form 

cl>{P^,x) = A,ix){Y,ix)r + A2{x)iY2ix)r 

and 

0(a,x) = (Fi(x)r + (F2(x)r-2. 

Proof. By Lemma 2.4, we notice that 0(P^,x) satisfy the recursive formula f{n,x) = 
xf{n — — f{n — 2,x). Therefore, the general solution of this linear homogeneous 
recurrence relation is f{n,x) = Di{x){Yi[x))"' + D2{x){Y2{x))"' . By some elementary 
calculations, we can easily obtain that Di{x) = Ai{x) for (f){P^,x), i = 1,2, from the 
initial values (j){P^,x), (f){Pj,x). 

By Lemma 2.4, (piCn, x) = (p{Pn, x) — (p{Pn-2, x) —2 and 4>{Pn, x) satisfy the recursive 
formula f{n, x) = f{n — l,x) — f{n — 2, x)—3. Similarly, we can obtain the general solution 
of this linear nonhomogeneous recurrence relation from the initial values (j){Pi,x) = x, 
0(P2,x) = a;2 - 1. □ 

Proof of Theorem 1.3 For n = 8, 12, 14, it is easy to verify E{P^) > E{Cn)- In 
the following, we always suppose n > 16. Using Lemma 2.2, we can deduce 



EiCr.) - E{P^ 



1 



TT 



log 



4>{Cn,ix) 



(P{Plix) 



dx. 



From Lemma 2.5, we have 

ix) = {Y,{ix)r + {Y2{ix)r - 2 = {{Z2{x))\x^ + 1) - {Z2{x)fx){Z,{x)r ■ 
+ {{Z^{x)Y{x' + 1) - {Z^{x)fx){Z2{x)r ■ ^" - 2, 
0(P„^ ix) = Ai{ix){Y^{ix)r + A2{ix){Y2{ix)T 



'Mx)f + zl 



{Z2{x)f + {Z2{x)y^ 



Firstly, we will prove that E{Cn) —E{P^) is an decreasing function of n for = 4k+j, 
j = 1,2,3, namely. 



log 



Ai{tx)iY,{tx))-+^ + A2itx)iY2itx)) 



n+4 



log 



Ai{ix){Yi{ix))'' + v42(zx)(F2M)" 



log 1 



Ko{n,x) 
Ho{n,x) 



< 0. 



A 



Case 1 n = 4:k + 2. 

In this case, Ho{n,x) = \(f){Cn,ix) ■ </)(P^_,_4, za:)| > and 

Ko{n,x) = {A,{ix) - A2(zx)) ((^2^)^ - (Fi(zx))^) - 2Ai(za;)(Fi(za;))"(l - {Y^{ix)Y) 
-'lA-,{ix){Y2{ix)r{\-{Y2{ix)f). 
Then, by some elementary calculations, we have 

K^in, x) = + l)(x^ + + 30x^ + 46x^ + 28x 

+ (Z2(x))"(x^ + 5x^ + 6x + Vx^ + A {x^ + 3x^ + 4)) 
+ + + 6x - Vx2 + 4 {x^ + 3x^ + 4))). 

If X > 0, then Zi{x) > 1, — 1 < Z2{x) < 0, and we obtain 

Ko{n,x) = x{x? + l)(Zi(a;))"g(n,a;) < + l)(Zi(a;))"g(10, x), 

where 

g(n,x) =(Z2(x))"(x^ + + 30x^ + 46x^ + 28x) 



+ (Z2(x))2"(x^ + Sx^ + 6x + Vx2 + 4 (x^ + ^ 4^,^, 
+ x^ + 5x^ + 6x - Vx2 + 4 (x" + 3x2 _^ 
By some simplifications, 

g(10,x) = -^x(x2 + 4)(2x^ + 17x^ + 47x^ + 46x2 + 10)- 



X 



10 



+ lOx® + 35x^ + 50x^ + 25x2 + 2 - Vx^ + 4(x^ + 8x^ + 21x^ + 20x^ + 5x)). 
Since 

(x^° + lOx^ + 35x^ + 50x^ + 25x2 + 2) ^ - (^Vx2 + 4 (x^ + 8x^ + 21x^ + 20x^ + 5x)) ^ = 4, 

we have g(10, x) < 0, and hence < 0. Similarly, we can prove < for x < 0. 

Therefore, we have shown that E{Cn) — E{P^) is an decreasing function of n for 
n = 4A; + 2. 

Case 2n = Ak + j,j = 1, 3. 

In this case, Ho{n,x) = |0(C„,zx) ■ 0(P^_,_4, zx) | > and 

Ko{n,x) = |((Fi(zx))"+^ + {Y2{tx)r+^ - 2)(Ai(zx)(Fi(zx))" + A2{ix){YSx)T)\ 

- |(Ai(zx)(ri(^x))"+^ + A2{ix)iY2{ix)r^'){iYSx)r + iY2{ix)r - 2)1 



= A/p(n,x) - v^w(n,x), 

where 

p(n,x) = (A2(«x)(Zi(x))^ + Ai(2x)(Z2(x))^ - Ai(2x)(Zi(x))2"+^ - A2(2x)(Z2(x))2"+4)' 
+ (-2Ai(zx)(Zi(x))" - 2A2(2X)(Z2(X))")' , 

w{n,x) = [ASx){Z^{x)Y + ASx){Z2{x)Y - A,{ix){Z^{x)f'-+^ - ASx){Z2{x)f''+'f 
+ {-2A,{ix){Z,{x)r+^ - 2A2{ix){Z2{x)r^^Y . 



Now we only need to check p{n,x) — w{n,x) < for all x and n. Firstly we suppose 
n = Ak + l. If X > 0, then > (^2(a;))2'^ < (^2(a;))^°, and we have 

p(n, x) - w{n, x) = x[x^ + 2)^(x^ + \f[x}^ + llx^ + 46x^ + 92x^ + 88x^ + 28x 



- 2(Zi(x))2"(v/;^^2^ (3,2 + 2) + x) 
+ 2(Z2(x))2"(v/^^(x2 + 2) - x)) 

< p(5, x) — w(5, x) 

= - x^{x^ + 4)(x2 + lf{x^ + 2)3(2x^ + 19x^ + 60x^ + 68x^ + 14) < 0. 

If X < 0, then (Zi(x))2" < {Z^{x)f^ , (^2(3;))^" > (^2(3;))^°. Similarly, x) - w(n, x) < 
p(5, x) — ^^(S, x) < 0. By the same discussion as the case of n = 4A; + 1, for = 4/c + 3 
and X > or X < 0, we can deduce that 

•p{n, x) — w{n, x) < p(7, x) — u;(7, x) 
= -x2(x2 + 4)(x2 + 2)3(x2 + l)3. 

(2x^'^ + 30x^2 + 178x^° + 533x® + 849x^ + 690x^ + 242x2 + 22) < 0. 

Thus, we have done for n = 4/c + j, j = 1, 3. 

Therefore, we have shown that E{Cn) — E{P^) is an decreasing function of n for 
n = Ak+j, J = 1,2,3. So, when n = 4A; + 2, E{Cn) - E{P^) < E{Cis) - E{P^s) = 
-0.03752 < 0; when n = Ak + 1, E{Cn) - E{P^) < E{Cn) - E{Pfj) = -0.00961 < 0; 
when n = 4:k + 3, E{Cn) - E{P^) < E{Ci^) - E{Pfg) = -0.02290 < 0. 

Finally, we will deal with the case of n = Ak. Notice that in this case both (f){Cn, ix) 
and (piPn, ix) are polynomials of x with all real coefficients. When n — j- 00, 

(Fi(ix))" + (F2(ia;))" - 2 ^ / Alfc) if a; > 



Ai{tx){Yi{tx))^ + A2{tx){Y2{tx))^ 1 ifx<0. 

' A2(tx) 



In this case, we will prove 



(ri(zx))"+(r2(^a;))"-2 , 1 

log — — -, TT ; r— < log ■ 



Ai(?x)(Fi(ix))" + A2{ix){Y2{ix))'' Ai{ix) 
for X > and 

(ri(ix))" + (^2(^0;))" -2 , 1 



Ai{ix){Yi{ix))"' + A2{ix){Y2{ix))'' A2{ix) 
for X < 0. In the following we only check the case of x > as the case of x < is similarly. 
Assume 

, iY,{ix)r + {Y2itx)r-2 1 / K,{n,x) 

log . N/.^/- NN^ , . N/.^/- NN. - log -TT-T = log 1 + 



Ai{ix){Yi{ix))'' + A2{ix){Y2{ix))'' Ai{ix) \ Hi{n,x) 

by some elementary calculations, we obtain Hi (n, x) > and 

[x^ + 9x^ + 28x^ + 36x2 + 16 + ((Z2(x))'" - l)Vx2 + 4(x^ + 7x^ + 16x^ + 14x)' 

< -^V^ ■ (x^ + + 28x^ + 36x2 + ig _ Vx2 + 4 (x^ + 7x^ + 16x^ + Ux)] < 0, 
x2 + 4 V V 



ct 



since 



{x^ + 9x^ + 28x^ + 36x2 + 16)^ - (Vx^ + 4 (x^ + + IGx^ + Mx)' 
4x^ + 48x^ + 204x^ + 368x^ + 256 > 0. 



Notice that if x > 0, then Ai{ix) = ^Zi{x)Y^+z'' ^ if x < 0, then A2{ix 

^^(")^^+^^ - ^i(^)-(^2(x)/8+/7) _ h-z^(x)h ^ 0. Thus, by Lemma 2.1, we have 



{Z2{x)f+Zl Z^(xy{{Z2(x)f+Zl) {Z2{x)r+z^ 

1 /'+°° 1 1 /■+°° / 1 

log — -— dx < - / -TT-T - 1 1 dx = -0.047643; 



1 /"^ 1 1 /'+°° / 1 

log , , , dx < - / , , , - 1 1 dx = -0.047643. 



Therefore, 



1 Z""*""" 1 f 1 

EiCn) - E(P^) < - log , , , dx + / log , , , dx < -0.047643 - 0.047643 < 0. 

TT Jo v4i(ix) A2(zx) 

The proof is completed. □ 
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